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一 种 求解 动态 优化 问题 的 免疫 文化 基因 算法 
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摘 要 : 针对 传统 免疫 网 络 动态 优化 算法 局 部 寻 优 能 力 弱 、 寻 优 精 度 低 及 易 早 熟 收敛 的 缺点 ， 提 出 一 种 求解 动态 优化 
问题 的 免疫 文化 基因 算法 。 基 于 文化 基因 算法 基本 框架 ， 将 人 工 免 疫 网 络 算 法 作为 全 局 搜索 算法 ,采用 禁 忆 搜索 算法 
作为 局 部 搜索 算 子 ; 同时 引入 柯 西 变异 加 强 算法 的 全 局 搜索 能 力 ， 并 有 效 防止 早熟 收敛 。 通 过 对 经 典 动 态 优 化 函数 测 
试 集 在 相同 条 件 下 的 实验 表明 ， 该 免疫 文化 基因 算法 相 较 于 其 他 同类 算法 具有 和 较 好 的 搜索 精度 和 收敛 速度 。 
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Immune-based memetic algorithm for dynamic optimization problems 


Yang Zhoul, Yuan Yichuanl, Luo Tingxing’, Qin Jin! 
(1. College of Computer Science & Technology Guizhou University, Guiyang 550025, China; 2. Guiyang Information Industry 
Development Center, Guiyang 550081, China) 


Abstract: The traditional immune network optimization algorithm has the shortcomings of weak local searching ability, low 
precision and premature convergence. In order to improve the algorithm performance, this paper proposed an artificial-immune- 
network-based memetic algorithm for dynamic optimization problems. Based on the framework of memetic algorithm, an 
artificial immune network algorithm serves as the global search algorithm, and a tabu search algorithm serve as the local search 


operator. At the same time, the algorithm introduced the Cauchy variation to improve global searching ability and prevent 


premature convergence. The experimental results show that the improved algorithm has better search precision and convergence 
speed compared with other algorithms. 
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8 引言 por=|® XxeX(tf csS,teT 如 
在 现实 世界 中 ， 许 多 优化 问题 往往 具有 动态 变化 的 性 质 。 5 表示 搜索 空间 ，1 表示 时 间 ， 了 表示 目标 函数 ，x 为 可 
例如 调度 问题 中 ， 新 任务 的 不 断 到 达 ， 产 品 原材料 质量 的 不 断 ” 行 解 ，X (GD) 是 在 1 时 间 的 可 行 解 集合 ，f (x,) 表示 在 1 时 刻 候 
变化 ,机 器 随时 间 增 加 的 不 断 磨损 以 及 随时 可 能 发 生 的 故障 等 。 ” 选 解 的 目标 函数 值 。 

也 就 是 说 ， 动 态 优化 问题 中 某 些 目标 或 条 件 随 着 时 间 不 断 发 生 动态 优化 问题 与 静态 优化 问题 的 区 别 在 于 静态 优化 问题 只 
着 变化 ， 某 个 时 刻 最 坏 的 选择 在 下 一 时 刻 有 可 能 会 变 好 ， 而 有 需要 算法 收敛 到 全 局 最 优 值 即 可 ， 而 动态 优化 问题 不 仅 需要 算 
些 条 件 本 来 无 关 紧 要 却 突 然 变 得 重要 。 社 会 变化 日 新 月 异 ， 动 ”法 找到 每 个 时 刻 的 全 局 最 优 值 ， 而 且 还 需要 具备 快速 追踪 随时 
态 优化 问题 所 涉及 领域 已 经 扩大 到 包含 几乎 所 有 工程 、 科 学 、 变化 的 最 优 解 的 能 力 ， 这 就 要 求 动 态 优化 算法 兼 具 快速 收敛 和 
经 济 和 生活 方面 的 问题 。 因 此 ， 对 于 动态 优化 问题 及 其 解决 方 保持 多 样 性 的 能 力 *** 玫 目前, 解决 动态 优化 问题 的 算法 
法 的 研究 越 来 越 受到 研究 人 员 的 关注 有。 主要 包括 进化 算法 和 群 智 能 算法 两 类 。 进 化 算法 是 基于 生物 进 
在 目前 的 研究 中 ， 动 态 优化 问题 通常 被 视 做 一 系列 静态 优 ”化 机 制 的 一 类 优化 算法 。 传 统 进化 算法 通常 用 于 求解 静态 优化 
化 问题 的 组 合 ， 其 动态 性 主要 表现 在 环境 变化 的 频率 以 及 环境 问题， 研究 者 通过 加 入 多 样 性 、 记 忆 机 制 、 预 测 机 制 、 多 种 群 
变化 的 剧烈 程度 。 动 态 优化 问题 9 可 以 定义 如 下 : 等 策略 ， 将 其 改进 为 求解 动态 优化 问题 的 进化 算法 。 群 智能 人 
法 是 模仿 生物 群体 的 行为 和 相互 间作 用 的 一 类 优化 算法 ， 包 括 
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蚁 群 优化 算法 


优化 算法 铺 识 ! 示 找到 引用 源 ，、 
其 衍生 出 一 ER 
就 是 受到 烟花 爆炸 的 启 
物 免 疫 系统 功能 的 一 类 智能 方法 ,广泛 应 用 于 
控制 等 领域 。 人 工 免疫 算法 的 记忆 机 制 、 


人 AAA AGRISAA 


上 错误 ! 未 找到 引用 源 。、 


Tt 


二 


ee 
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疫 算法 可 改进 范围 
结合 。 例 如, Liu 等 
工 免 疫 网 络 算法 (CoAIN) 和 RAR 
作 性 ， 使 网 络 中 个 体 不 仅 相互 
算法 搜索 能 力 。 赵 辉 等 人 提出 一 种 
人 工蜂 群 免疫 算法 宽 
算法 , 克服 传统 人 工 由 


a 
制 ， 而 且 相 互 合作 ， 从 了 
6 式 精英 进化 模型 的 
蜂 群 思想 和 免疫 克隆 选择 
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， 继 续 ， 否 则 转 入 h); 
缚 细胞 抑制 ,抑制 距离 小 于 设 定 值 的 两 细胞 个 体 中 适应 度 低 
细胞 个 体 ; 
多 引 入 新 细胞 : 按 比例 随机 产生 新 细胞 个 体 加 入 当前 细胞 


h) 终 止 条 件 判断 :判断 当前 是 否 达到 终止 条 件 ， 若 达到 终 
条 件 ， 输 出 结果 ， 否 则 ， 转 入 b)。 
1.2 ”禁忌 搜索 算法 

禁忌 搜索 算法 (tabu search,TS) 是 一 种 元 启发 式 搜索 算法 
“， 通过 模仿 人 类 的 记忆 功能 , 使 用 禁忌 表 来 封锁 刚刚 搜 
区 域 来 避免 迁 回 搜索 ， 从 而 使 算法 能 够 跳出 局 部 极 值 。 
也 传统 优化 算法 ， 禁 忌 搜 索 算 法 具有 更 强 的 “ 拒 山 ”能 
部 搜索 能 力 。 禁忌 搜索 算法 在 旅行 商 问 题 、0-1 背包 问题 
合 优化 问题 上 已 有 十 分 广泛 的 应 用 ， 近 年 来 许多 学 者 也 演 


运用 于 函数 优化 问题 的 求解 。 


文化 基因 


算法 (memetic algorithm,MA) 是 近 些 年 演 


领域 的 一 个 研究 热点 ， 其 将 全 局 搜索 与 局 域 搜索 结合 的 思想 有 
利于 平衡 动态 优 与 开发 间 的 关系 。 因 此 ， 本 文 将 在 
已 有 研究 的 基础 之 上 ， 提 出 一 种 免疫 文 
at 
en 
在 算法 寻 优 到 一 
前 最 优 解 邻 域内 进行 搜索 ， 提 高 


精度 低 、 


i 


禁忌 搜索 算法 其 本 质 是 一 种 优秀 的 邻 域 搜索 算法 ， 具 有 很 
快 的 收敛 速度 。 但 其 收敛 速度 的 快慢 和 寻 优 结果 的 优 劣 很 大 程 
度 上 取决 于 初始 解 的 好 坏 ， 好 的 初始 解 往往 能 够 帮助 算法 快速 
上 全 局 最 优 解 ， 而 差 的 初始 解 将 可 能 大 大 降低 算法 收敛 速 


。 禁 忌 搜 索 算法 的 关键 在 于 邻 域 的 选取 、 禁 忌 表 的 设计 以 及 


邻 域 体现 局 部 搜索 思想 ， 禁 忌 表 体 现 算法 避免 迁 回 


特点， 特赦 准则 是 对 优秀 解 的 一 种 反馈 。 禁 忌 搜 索 算法 


因 陷 入 局 部 最 优 而 售 


。 同 时 ， 引 入 柯 西 变异 算 村 


变异 算 子 切换 使 用 ， 


实验 仿真 验证 


1 ”相关 研究 


1.1 dopt-aiNet 算法 
人 工 免 疫 网 络 算法 


法 全 局 寻 优 能 力 。 通 过 对 函数 集 世 
所 提出 的 免疫 文化 基因 动态 优化 算法 的 性 能 。 


于 生物 免疫 系统 的 信息 处 


起 来 ， 主 要 涉及 
疫 学 原理 。 A gil 
缩 和 数据 聚 类 。 
法 求解 优化 问题 ， 
3 之 后 ， 又 在 opt-aiNet 的 基础 
网 络 动 态 优化 算法 (dopt-aiNet) 生 相交 到 引用 滨 , 。 并 在 
化 竞赛 ， 在 求解 动态 优 


CEC 举办 的 动态 优 


性 台 已 错误 ! 未 找到 引用 源 。 


dopt-aiNet 算法 模仿 生物 免疫 细胞 发 现 外 来 抗原 ，i 


Fe 


六 过 


et 
提出 人 工 免 疫 网 络 优 


化 算法 (opt-aiNet) 生 于 丰 并 到 优 解 了”Y， 继 续 ， 否 则 直接 转 入 6); 
已 


思路 主要 是 搜索 当前 一 个 解 的 邻 域 ， 如 果 解 邻 域内 的 最 
tC 解 优 于 当前 解 ， 将 当前 解 替换 为 该 解 邻 域内 的 最 优 解 ， 并 继 
在 当前 解 的 邻 域 内 进行 搜索 , 直到 算法 终止 。TS 算法 流程 可 


几 个 步骤 : 


a) 初 始 化 参数 :给 定 初始 解 和 ,计算 目标 函数 值 f(X,) ， 
将 初始 解 赋值 给 当前 解 克 =X。; 

b) 生 成 邻 域 和 候选 解 集 : 生成 当前 解 的 邻 域 N(X)， 从 
中 选 出 一 定数 量 的 解 作为 候选 解 集合 ， 并 计算 所 有 候选 解 的 


满足 邻 域 搜索 条 件 : 如 果 候 选 解 集合 不 空 并 且 


忌 准 则 : 从 候选 解 集合 中 找到 最 优 解 了 ， 如 果 了 了 


在 禁忌 表 中 ， 将 了 赋值 给 当前 最 优 解 ， 否 则 ， 将 其 从 当前 解 
删除 ， 转 入 9); 
禁忌 表 更 新 : 如 果 找 到 最 优 解 了 了 ， 将 当前 解 X 插入 禁忌 
武 值 X = 了 了; 


原 ， 刺 激 免 疫 


消灭 抗原 的 免疫 过 程 。 


计算 细胞 个 体 适应 度 ; 


b) 细 胞 克隆 ,对 每 个 细胞 


dj) 精英 选择 ,在 每 组 变异 后 的 细胞 


其 具体 流程 如 1 
a) 初 始 化 种 群 , 生 
c) 细 胞 突变 ,对 克隆 

高 的 细胞 个 体 留 下 ; 


9) 细 胞 抑制 


有 判断 结束 : 如 果 满 足 结束 条 件 , 则 终止 计算 , 否则 转 入 b)。 
免疫 文化 基因 动态 优化 算法 


基因 算法 是 基于 Darwinian 的 自然 进化 理论 和 Darkins 


进化 思想 产生 的 一 种 混合 算法 框架 下。 算法 将 
进化 算法 和 某 些 局 域 搜索 方法 结合 在 一 起 ， 进 化 算法 负责 执行 


判断 ,判断 细胞 种 群 是 否 达到 抑制 条 件 


全 局 搜索 ， 局 部 搜索 策略 用 于 改善 某 些 个 体 的 搜索 精度 ， 从 而 
法 在 探索 和 开发 间 保持 习 


I 


衡 。 
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2.1 算法 思想 

人 工 免疫 网 络 算法 (dopt-aiNet) 在 求解 动态 优化 问题 时 ， 主 
要 通过 初始 化 、 克 隆 扩张 、 高 频 变 异 、 精 英 选 择 等 操作 来 优化 
候选 解 ， 最 终 找 到 问题 的 解 ， 并 且 通 过 抑制 机 制 和 抑制 之 后 引 
入 新 个 体 的 方式 维持 解决 动态 优化 问题 所 必须 的 种 群 多 样 性 。 
算法 种 群 个 体 广 泛 分 布 在 解 空间 各 个 区 域 往往 能 体现 良好 的 全 
局 搜索 能 力 ， 但 当 算法 寻 优 到 一 定 程度 时 ， 容 易 出 现 寻 优 能 
不 足 ， 导 致 早熟 收敛 ， 一 定 程度 上 影响 了 算法 的 寻 优 精度 。 针 
对 dopt-aiNet 算法 求解 动态 优化 问题 存在 的 不 足 ， 本 文 基于 文 
化 基因 算法 思想 提出 解决 动态 优化 问题 的 免疫 文化 基因 算法 
(artificial immune network based memetic algorithm, AINMA)。 利 
禁忌 搜 索 的 局 部 搜索 优势 ， 结 合 人 工 免疫 网 络 算法 全 局 搜索 
E 力 强 的 特点 ， 通 过 在 细胞 抑制 之 后 引入 禁忌 搜索 ， 人 工 免疫 
络 为 禁忌 搜索 提供 好 的 初始 解 , 提高 算法 的 搜索 精度 ,此 外 ， 
入 柯 西 变异 算 子 与 传统 高 斯 变异 算 子 结合 变异 ， 寻 优 前 期 使 
j 柯 西 变异 保证 更 大 限度 的 搜索 整个 解 空间 ， 寻 优 后 期 使 用 高 
斯 变异 算 子 保证 算法 搜索 精度 。 同 时 ， 改 进 算 法 采用 设立 侦 测 
细胞 的 方式 侦 测 环境 变化 ， 侦 测 种群 在 每 次 迭代 之 后 计算 适应 
度 ， 如 果 适 应 度 发 生变 化 即 环境 发 生 改变 ， 引 导 寻 优 种 群 重新 
开始 搜寻 新 的 最 优 值 。 
2.2 柯 西 变异 算 子 

传统 高 斯 变异 算 子 是 对 个 体 ,= (2 ) 按照 式 (2) 
进行 操作 。 


zw 


习 


3 du 


X=% +t7*G(0,D O) 
其 中 : 7 为 控制 变异 步 长 的 常数 ，C(0,1D 是 均值 为 0, 偏差 
为 1 的 高 斯 分 布 函数 产生 的 随机 数 。 
柯 西 变 异 算 子 是 对 个 体 X, = Co ) 按照 公式 (3) 
进行 操作 。 


x =X + *C(0,1) (G3) 

其 中 :7 为 控制 变异 步 长 的 常数 ，C(0,D 是 由 t=1 的 柯 西 

分 布 函 数 产生 的 随机 数 。 
比较 高 斯 密度 函数 和 柯 西 密度 函数 可 知 ， 柯 西 密度 函数 两 

端 尾 巴 比 高 斯 密度 函数 更 长 ， 这 就 使 得 由 柯 西 变异 产生 的 变异 
子 代 与 父 代 差异 更 大 ,个 体 则 有 更 大 几率 跳出 局 部 最 优 , 因此 ， 
柯 西 变异 往往 比 高 斯 变异 具有 更 好 的 全 局 搜索 能 力 。 本 文 将 两 


杨 洲 ， 


合作 期 下 
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放 最 近 访 问 过 的 邻 域 候选 解 ， 从 而 防止 迁 回 搜索 ， 也 可 以 避免 


算法 陷入 局 部 最 优 。 


禁忌 搜索 中 特赦 ; 


AINMA 算法 中 , 当 某 个 邻 域 候选 解 的 适应 度 优 于 
无 论 这 个 候选 解 是 否 在 禁忌 表 中 ， 都 接受 这 个 候选 解 。 


全 则 的 设立 是 为 了 防止 错过 优良 解 。 在 


历史 最 L 解 ， 


禁忌 搜索 中 , 邻 域 的 划分 十 分 关键 ,不 同 于 组 合 优化 问题 ， 
函数 优化 问题 的 解 空 间 是 连续 的 ， 因 此 ， 本 文采 用 文献 错误 ! 未 

找到 引用 源 。 中 提 及 的 邻 域 划分 方式 选取 邻 域 解 。 
将 动态 优化 问题 中 单个 静态 优化 问题 用 如 式 (4) 的 模型 表 


minf(X) st. L<X<U (4) 


其 中 : 和 为 n 维 变量 ，XX = (X64, 鸭 ,…, 鸭 ) ， 对 应 分 变量 的 


上 下 限 民 分 别 为 上 = (l,l 2， “sb 3 U =(U,U,,,,) o 


同心 超 和 矩阵 将 当前 解 的 分 量 如 (7 二 上 …,1) 分 成 个 空 


间 。 每 个 空间 按 式 (5) 划 分 。 


R(x,r 117)= {x | | 


lL< x < j= 2 


n= *2,1=1,2,.…,k (5) 


R(x,r)={x'"| Ey —X, 


1<X <u,j=1,2,.…,n 
4 J 


在 RR 的 每 个 同心 矩形 内 随机 选取 一 个 点 , 选 出 的 上 个 点 组 
成 当前 解 分 量 x 的 邻 域 候选 解 {Xi, 台 ,…, 和 }。 


2.4 算法 流程 


AINMA 算法 可 分 为 以 下 几 个 步骤 : 
a) 初 始 化 种 群生 成 寻 优 种 群 和 使 测 种 群 ， 并 设置 相应 参 


数值 ; 


b) 克 隆 : 寻 优 种 群 
0c) 变异 : 对 每 个 克 


个 体 克 隆 扩张 ; 
隆 个 体 进行 高 频 突变 ， 如 果 变 异 标志 为 


0, 选用 柯 西 变异 算 


进行 突变 , 否则 选用 高 斯 变异 算 子 进行 突 


变 ， 并 计算 变异 个 体 的 适应 度 值 ; 


dd) 选择 :在 每 组 变异 个 体 中 选择 适应 度 大 的 个 体 留 下 ; 
6) 抑制 判断 : 判断 是 否 达 到 抑制 条 件 , 若 达 到 条 件 , 继续 ， 


否则 转 入 D; 


抑制 : 变异 标志 设置 为 1， 抑制 距离 小 于 设 定 值 的 两 个 体 


种 变异 方式 结合 使 用 ， 在 算法 寻 优 前 期 ， 使 用 柯 西 变异 能 够 防 
止 寻 优 个 体 陷入 局 部 最 优 导 致 算法 提前 停滞 ， 当 算法 满足 抑制 
条 件 进入 抑制 过 程 ， 也 就 是 说 寻 优 已 经 到 达 一 定 程度 ， 可 能 十 
分 靠近 全 局 最 优 ， 此 时 使 用 高 斯 变异 防止 过 度 变 异 导致 寻 优 个 
体 远离 全 局 最 优 值 。 
2.3 ”禁忌 表 、 邻 域 与 特赦 准则 
本 文 将 禁忌 搜索 融合 于 dopt-aiNet 算法 中 ， 利 用 免疫 算法 
克隆 、 变 异 、 选 择 等 操作 为 禁忌 搜索 提供 一 个 较 好 的 初始 解 ， 
再 通过 禁忌 搜索 作为 局 部 搜索 工具 提高 整个 算法 的 寻 优 精度 。 
以 下 是 对 AINMA 算法 中 禁忌 表 、 邻 域 和 特赦 准则 的 简要 介绍 : 
禁忌 表 是 禁忌 搜索 算法 的 核心 ， 在 AINMA 算法 中 用 来 存 


中 适应 度 低 的 个 体 ; 


多 禁忌 搜索 : 生成 


Ik 


4 前 最 优 解 的 邻 域 和 候选 解 集合 ， 在 候 


选 解 集中 寻找 最 优 解 ， 并 不 断 更 新 禁忌 表 ， 重 复 该 搜索 过 程 ] 
到 满足 结束 条 件 ; 


h) 引 入 随机 个 体 : 


四 


随机 产生 新 个 体 加 入 当前 寻 优 种 群 中 ; 


iD 侦 测 环境 变化 : 利用 侦 测 种 群 侦 测 环境 是 否 发 生变 化 , 若 
发 生变 化 ， 将 变异 标志 设置 为 0， 和 否则， 继续 ; 


停止 判断 : 判断 当前 是 否 达到 停止 条 件 , 若 达到 停止 条 件 ， 
输出 结果 ， 否 则 ， 转 入 b)。 


AINMA 算法 流程 


图 如 图 1 所 示 。 
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Ky 人 化 基因 


录用 稿 杨 洲 ， 
(= 和 王公。 mi IA 
me 3 ”实验 与 分 析 
| 初始 化 种 群 ， 计 算 种 群 适应 度 | 
ET 3.1 实验 设置 
OO 克隆 扩张 有 、 二 、 > 各 
为 了 验证 本 文 提出 的 免疫 文化 基因 动态 优化 算法 的 性 能 ， 
是 一 一 一 否 
何 西 变异 ] 而 出 彼 萤 | 采用 文献 错误 ! 未 找到 引用 源 。 中 提出 的 generalized dynamic 
计算 种 凌 天 应 度 benchmark generator(GDBG) 浮 数 测试 集 进行 仿真 测试 。GDBG 
函数 测试 集 包含 6 种 测试 函数 类 型 ， 每 种 测试 函数 对 应 7 种 改 
变 类 型 ， 包 括 小 幅 改变 、 大 幅 改变 、 随 机 改变 、 混 沌 改变 、 周 
下 
asial 期 性 改变 、 周 期 带 噪声 改变 和 维度 改变 。 该 动态 函数 测试 集 | 
人 表 1 所 示 5 种 函数 通过 旋转 、 组 合 形成 以 下 6 类 不 同 难度 的 问 
[生成 当前 最 优 解 邻 域 和 候选 解 集 
各 了 题 : 
下 | 在 候选 解 集中 寻找 最 优 解 | 
六 二 Fl: rotation peak function(with 10 and 50 peaks) 
os F2: composition of Sphere’s function 
F3: composition of Rastrigin’s function 
一 F4: composition of Griewank’s function 
FS$: composition of Ackley’s function 
| F6: hybrid composition function 
Fl 是 最 大 值 优化 问题 ，F2-F6 是 最 小 值 优化 问题 。 其 中 ， 
图 1 AINMA 算法 流程 图 Fl 分 为 两 种 测试 函数 ， 一 种 带 有 10 个 峰 ， 另 一 种 带 有 50 个 
峰 
表 1 GDBG 基础 组 成 函数 表达 式 
name function range 
Sphere f= [-100,100] 
Rastrigin f= 0 —10cosQm,)+10) [-5,5] 
ee a 
XxX) = 和 [a”* cos(2np” x;+O0.S)MD—n [a”* cos(ap”)] 
Weierstrass A 之 - (名 9 和 [-0.5,0.5] 
a=0.5,b=3,k,,.. =20 
1 n 2 n Xx, 
1 XxX) 一 一 -一 一 (Xx,) 一 CoOS( 一 一 ) 十 1 a 
Griewank Oe > Co) 一 [ [cos' 7 [-100,100] 
Ackley fx) =—20exp(—0.2,|1 > x?)—exp(t DcosQm)) +20+e [-32,32] 
i il 
本 文选 择 使 用 误差 均值 (Average mean) 作 为 算法 性 能 评价 实验 中 ， 测 试 函数 的 参数 设置 包括 维度 d =10 ， 变 化 范 辕 


指标 。 计 算 方式 如 式 (6) 所 示 。 
mean = > > E™ (nD)/(runs*num_change) (6) 
其 中 : runs 表示 算法 独立 运行 的 次 数 ，num _change 表示 


每 种 改变 发 生 的 次 数 ，E” (1) 表示 算法 在 第 i 次 独立 运行 中 第 


j 次 改变 时 的 适应 度 绝对 误差 值 ， 如 式 (7) 所 示 。 
E™(D)=|f 60)) -fx 0) 0) 


其 中 : f(x (?)) 和 (x (D) 分 别 表示 算法 寻 得 最 优 值 和 
里 论 最 优 值 。 


d e[5,15] ; 搜索 范围 x e[-5,5] ; 改变 频率 fre =10000*4d ，; 
改变 次 数 num _ change = 60 ; 变化 周期 period =12 。 算法 的 
参数 设置 包括 : 种 群 规模 和 =10; 克隆 规模 C = 3 ;抑制 阔 值 
6 =5; 变异 步 长 7 =1。 
3.2 ”结果 与 分 析 

为 验证 本 文 提出 的 免疫 文化 基因 算法 (AINMA) 的 性 能 ， 利 
用 基本 人 工 免疫 网 络 动 态 优化 算法 (dopt-aiNet)、 粒 子 群 动态 优 
化 算法 (CPSO)、 变 量 重 定位 正 交 粒子 群 算法 (OLPSO-VRS)、 免 
疫 文化 基因 算法 (AINMA) 以 及 不 含 禁忌 搜索 的 免疫 文化 基因 
算法 (AINMA without TS) 对 GDBG 动态 函数 测试 集 进 行 了 相同 
条 件 下 的 仿真 实验 。 结果 对 每 个 测试 函数 进行 5 次 独立 测 


实验 


试 后 平均 得 出 。 仿 真 结果 如 表 2 所 示 ， 表 中 数 


居 均 以 科学 计数 


杨 洲 ， 等 : 一 种 求解 动态 优化 问题 的 免疫 文化 基因 算法 


表示 ,并 将 4 种 算法 在 每 种 情况 下 的 误差 均值 最 小 值 进 行 加 粗 。 
表 2 AINMA、AINMA(without TS)、dopt-aiNet、CPSO、OLPSO-VRS 对 GDBG 的 误差 均值 对 比 


Func. Algorithm Tl T2 T3 T4 TS T6 T7 
AINMA 2.76e-04 2.4Se+00 2.7Se+00 S.91e-04 1.91e+00 1.38e-01 2.29e+00 
AINMA(without TS) 5.81e-04 2.45e+00 2.97e+00 1.15e-03 2.26e+00 6.49e-01 3.17e+00 
ken dopt-aiNet 1.72e-01 4.57et+00 5.03e+00 6.34e+00 5.84e+00 1.03e+01 3.93e+00 
uy CPSO 4.78e-02 2.48e+00 4.69e+00 6.48e-02 1.39e+00 1.34e+00 3.78e+00 
OLPSO-VRS 2.08e-02 2.36e+00 3.57e+00 9.73e-02 1.36e+00 1.27e+00 3.54et+00 
AINMA 2.90e-04 1.47e+00 2.77e+00 1.07e-01 7.04e-01 2.86e-01 1.71e+00 
F1(50p AINMA(without TS) 5.76e-04 1.49e+00 3.59e+00 5.18e-01 7.75e-01 4.72e-01 2.31e+00 
dopt-aiNet 4.58e-01 4.41et+00 4.98e+00 3.09e+00 3.28e+00 5.94e+00 4.12e+00 
和 CPSO 3.74e-01 4.16e+00 4.84e+00 3.13e-01 2.01e+00 1.63e+00 3.77e+00 
OLPSO-VRS 1.91e-01 3.3Se+00 3.87e+00 1.27e-01 1.36e+00 9.34e-01 3.77e+00 
AINMA 2.94e-03 2.28e+01 2.21e+00 2.94e-03 3.46e+00 2.27e-01 7.76e-01 
AINMA(without TS) 5.76e-03 2.8Se+01 8.11e+00 2.15e-01 6.20e+00 3.27e-01 7.92e-01 
F2 dopt-aiNet 9.37e-01 1.08e+01 2.02e+01 1.86e+00 9.87e+01 5.99e+00 4.19e+00 
CPSO 1.17e+00 1.24e+01 8.54et+00 6.05e-01 2.37e+01 1.69e+00 4.12e+00 
OLPSO-VRS 8.82e-01 1.08e+01 8.28e+00 6.69e-01 1.66e+01 1.10e+00 2.03e+00 
AINMA 1.3Se+01 8.43e+02 8.47e+02 4.53e+02 8.01e+02 6.04e+02 5.26e+02 
AINMA(without TS) “2.16e+01 8.58et+02 8.91e+02 4.54e+02 8.02e+02 6.51e+02 5.26e+02 
F3 dopt-aiNet 7.65e+02 1.01e+03 1.00e+03 1.10e+03 1.04e+03 1.29e+03 1.11e+03 
CPSO 1.25e+02 8.04e+02 7.8Se+02 4.56e+02 8.36e+02 7.70e+02 5.97e+02 
OLPSO-VRS 1.5let+01 7.Sle+02 S.66e+02 3.22e+02 S.21e+02 3.79+02 4.1Se+02 
AINMA 3.06e-03 1.40e+01 6.78e+00 3.81e-03 .06e+01 6.71e-01 7.86e-01 
AINMA(without TS) 6.01e-03 5.48e+01 8.67e+00 6.57e-03 5.53e+01 6.84e-01 7.93e-01 
F4 dopt-aiNet 1.09e+00 1.3Se+02 8.68e+01 5.77e+00 3.01e+02 1.59e+01 5.59e+01 
CPSO 2.87e+00 3.47e+01 3.26e+01 7.85e-01 6.1Se+01 5.28e+00 1.08e+01 
OLPSO-VRS 1.24e+00 3.00e+01 2.18e+01 9.55e-01 4.13e+01 1.4Se+00 7.8Se+00 
AINMA 8.66e-03 1.61e-02 1.SSe-01 7.37e-03 2.54e+00 6.25e-01 1.80e-01 
AINMA(without TS) 1.84e-02 9.47e-01 1.55e-01 2.18e+00 3.79e+01 6.69e-01 2.02e-01 
FS dopt-aiNet 3.58et+01 3.67e+01 3.09e+01 1.17e+02 1.01et+03 4.57e+02 2.3Se+02 
CPSO 2.16e+00 2.65et+00 4.04e+00 1.12e+00 8.16e+00 .74e+00 5.89e+00 
OLPSO-VRS 5.64e-01 3.22e-01 4.52e-01 1.21e-01 S.99e-01 3.28e-01 6.20e-01 
AINMA 2.12e-01 2.34e+01 7.68e+00 2.97e+00 1.40e+01 4.40e+00 4.74e+00 
AINMA(without TS) 1.97e+00 5.04e+01 1.23e+01 4.30e+00 1.52e+01 4.53et+00 5.12et+00 
F6 dopt-aiNet 2.36e+01 3.48e+02 4.29e+02 8.67e+01 8.74e+02 5.27e+02 3.62e+02 
CPSO 7.93e+00 2.47e+01 2.86e+01 7.54e+00 6.79e+01 2.68e+01 3.48e+01 
OLPSO-VRS .72e+00 1.43e+01 1.25e+01 7.73e+00 1.94e+01 7.51let+00 9.77e+00 

根据 表 2 中 算法 对 比 结果 分 析 可 知 ， 本 文 所 提出 的 免疫 文 升 做 出 一 定 贡献 。 从 测试 函数 角度 分 析 , 由 于 F3 的 基础 组 成 函 
化 基因 动态 优化 算法 (AINMA) 在 所 有 49 种 测试 情形 下 ， 误 差 ” 数 Rastrigrin 函数 是 典型 非 线性 多 模 态 函数 ， 峰 与 峰之 间 起 伏 


均值 都 比 传 统 免 疫 网 络 动态 优化 算法 (dopt-aiNet) 小 。 也 就 是 说 ， 


较 大 , 并 且 具 


AINMA 算法 


索 的 算法 误差 更 小 ， 


有 更 强 的 寻 优 能 力 。 与 不 含 禁忌 搜索 的 免疫 文 
化 基因 算法 比较 , 在 大 多 数 情况 下 , AINMA 算法 比 不 含 禁忌 搜 
也 就 是 说 ， 禁 忌 搜 索 的 引入 为 算法 精度 提 


值 相 较 于 其 他 几 类 测试 函 
在 F3、F4、F4、F6 上 的 精度 提升 最 为 
杂 问 题 的 求解 能 力 有 很 大 改善 。 从 变 


大 量 


负极 值 点 ， 因 此 ,算法 在 F3 
数 大 。6 类 测试 函数 


上 的 误差 均 


hh ，AINMA 算法 
明显 , 表明 算法 对 于 较 复 
化 类 型 角度 分 析 ， 算 法 在 
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好 用 入 杨 洲 ， 等 : 一 种 rm YA 作 央 和 
Tl1、T4、T6 三 种 变化 类 型 下 的 误差 均值 相对 较 小 ， 即 算法 在 问 Wa [3 丰收 名 闭 国 
题 发 生 小 幅 变 化 、 混 沌 变化 和 周期 噪声 变化 情况 下 寻 优 能 力 较 
好 。 综合 分 析 , AINMA 算法 与 dopt-aiNet、CPSO、OLPSO-VRS | | 
三 个 算法 相 比 ，49 种 测试 情况 下 有 37 种 测试 情况 误差 均值 最 Ro | 
小 。 其中， 与 CPSO 相 比 ，AINMA 算法 仅 在 个 别 问题 上 不 如 本 | | | | ] 
CPSO 出 色 ; 与 OLPSO-VRS 相 比 ，AINMA 算法 仅 在 F3 上 稍 ! :A 
劣 于 OLPSO-VRS, 但 差距 并 不 明显 。 总 之 ，AINMA 算法 无 论 a 
在 收敛 精度 还 是 寻 优 能 力 上 都 有 极 大 改善 ， 并 表现 出 对 各 类 动 图 5 F3 在 Tl 收敛 曲线 图 
态 优 化 问题 广泛 的 适应 性 。 a TIE 
根据 仿真 实验 结果 ， 可 以 得 出 各 类 测试 函数 在 各 种 变化 类 。 
型 下 的 收敛 曲线 图 ， 以 下 分 别 列举 F1-F6 在 T1 变化 类 型 下 的 
收敛 曲线 图 进行 分 析 。 如 图 2 至 图 8 所 示 ， 各 图 纵 坐 标 表 示 适 起 | 
应 度 值 ,其 中 Fl 为 最 大 值 优化 问题 ,F2-F6 为 最 小 值 优化 问题 ， “| 1 | I! 
横 坐标 表示 环境 改变 的 次 数 ， 本 文 实验 设置 改变 次 数 为 60 次 ， wl | | i 1 
各 图 随机 选取 连续 10 次 变化 进行 分 析 。 一 一 一 一 一 一 
轩 PIG peaksy 人 ETI 履 从 让 图 
。 i | 面 面 图 6 F4 在 TI 收敛 曲线 图 
-Pp | | | | | Se 让 1] 收 化 作 急 同 
wb | | ee 
u 本 
人 i | | | | | 
图 2 F1(10 peaks) 在 TI 收敛 曲线 图 ”| | | | Eg 
Fl (G0 meaks) 在 T1 收 黎 卢 线 图 人 
100 gs | 改变 次 数 
20 i 和 二 图 7 F5 在 Tl 收敛 曲线 图 
加 民 | | ET 
， ， a " a | | | | | | | 
图 3 F1(50 peaks) 在 Tl 收敛 曲线 图 mj 1 || | ,| 
a ,克也 收 信 直 站 放 Pi : 人 5 
1000 改变 次 数 
四 图 8 F6 在 TI 收敛 曲线 图 
_ 在 动态 优化 问题 中 ,收敛 曲线 图 通常 反映 算法 的 收敛 速度 、 
wi | | | 1 对 于 环境 变化 的 响应 以 及 是 否 存在 提前 停滞 现象 。 通 过 观察 分 
中 2 一 二 析 算 法 在 各 类 问题 中 的 收敛 曲线 图 ， 一 方面 可 以 了 解 算法 对 该 
到 朗 交 娄 类 动态 优化 问题 的 适应 程度 ， 即 当 问 题 发 生变 化 时 ， 算 法 是 否 
图 4 F2 在 Tl 收敛 曲线 图 能 及 时 适应 环境 变化 ， 调 整 搜索 方向 ; 男 一 方面 可 以 了 解 算法 
在 每 次 变化 内 的 收敛 速度 ， 收 敛 速度 越 快 ， 算 法 越 能 够 适应 快 
速 变 化 的 环境 。 


从 图 2 可 以 看 出 ， 包 含 10 个 峰 的 Fl 测试 函数 在 Tl 变化 


尖 玫 


型 下 收敛 速度 较 1 


环境 发 生变 化 时 能 


岂 ， 收 敛 过 程 中 极 少 出 现 停 滞 状 态 ， 算 法 在 
及 时 响应 W 
所 示 是 包含 50 个 峰 的 Fl 测试 


a 图 3 
函数 在 T1 变化 类 型 下 的 收敛 曲 


于 峰 数 增加 导致 的 环境 复杂 性 变化 并 未 影响 算法 的 收 
图 4 表示 F2 测试 函数 在 Tl 
变化 类 型 下 的 收敛 性 ， 由 图 可 知 ， 算 法 第 一 次 从 较 差 初始 值 收 
不 境 发 生 改变 时 ， 算 法 的 适应 度 值 并 未 
大 幅 上 升 ， 表 明 算 法 始终 保持 良好 的 多 样 性 。 图 5 所 示 F3 测 
函数 本 身 的 复杂 性 ， 算 法 易 陷入 局 部 最 优 值 ， 并 且 
算法 对 于 环境 的 适应 速度 较 慢 ， 前 期 寻 优 结果 差 ， 但 当 算法 逐 
渐 适 应 环境 之 后 ， 其 寻 优 精度 不 断 提升 。 图 6 至 图 8 所 示 F4- 
F6 测试 函数 在 Tl 变化 类 型 下 的 收敛 曲线 图 所 反映 的 特征 与 图 
4 相似 。 
综合 以 上 分 析 可 知 ，AINMA 算法 继承 了 传统 人 工 免 疫 网 
络 算法 的 多 样 性 ， 并 在 此 基础 上 提高 的 收敛 速度 和 搜索 精度 ， 
当 环 境 发 生变 化 时 能 够 快速 响应 并 开始 新 的 搜索 。 分 析 表 明 ， 
AINMA 算法 对 于 求解 连续 函数 动态 优化 问题 有 良好 性 能 。 


mT 


4 ”结束 语 


本 文 在 传统 人 工 免疫 网 络 算法 (dopt-aiNet) 的 基础 之 上 ， 针 
对 算法 搜索 精度 低 和 寻 优 后 期 易 出 现 早 熟 收敛 的 缺陷 ， 提 出 了 
一 种 免疫 文化 基因 动态 优化 算法 (AINMA)。AINMA 算法 保留 
了 dopt-aiNet 算法 多 样 性 好 的 优点 ， 通 过 引入 柯 西 变 异 算 子 ， 
与 传统 高 斯 变异 算 子 交 蔡 使 用 ， 在 算法 不 同 寻 优 时 期 使 用 不 后 
的 变异 方式 ， 充 分 发 挥 两 种 变异 算 子 的 优势 ， 防 止 个 体 陷 入 局 
部 最 优 ， 一 定 程度 上 提高 了 算法 的 全 局 搜索 能 力 。 同 时 ， 利 1 
禁忌 搜索 的 局 部 寻 优 能 力 ， 将 其 融合 进 dopt-aiNet 算法 当中 ， 
有 效 提高 算法 的 寻 优 精度 ,通过 对 GDBG 动态 函数 测试 集 的 仿 
真实 验 表 明 ，AINMA 算法 相 较 于 改进 前 在 搜索 精度 和 收敛 速 
度 上 都 有 很 大 提高 ， 并 且 与 其 他 同类 动态 优化 算法 比较 也 有 明 


显 优势 。 

本 文 所 提出 的 AINMA 算法 对 于 存在 大 量 局 部 最 优 的 复杂 
动态 问题 表现 出 适应 性 差 的 现象 。 下 一 步 将 深入 分 析 ， 提 出 改 
进 。 
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